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EQUILIBRIUM NORMAL-SHOCK AND STAGNATION-POINT PROPERTIES 
OF HELIUM FOR INCIDENT-SHOCK MACH NUMBERS FROM 1 TO 30 

By Walter B. Olstad, Jane T. Kemper, and Roger D. Bengtson 
Langley Research Center 

SUMMARY 

Charts are presented of the thermodynamic properties of helium, including the 
effects of ionization and particle interactions, for temperatures to 100 000° K and pres- 
sures from 10-7 to 10^ atmospheres. The properties behind the incident, standing, and 
reflected shocks and at the in-flight and shock-tube stagnation conditions are also pre- 
sented for incident-shock Mach numbers from 1 to 30 and quiescent gas pressure from 
10"7 to 1 atmosphere. 


INTRODUCTION 

Helium has been used extensively in research facilities such as hypersonic wind 
tunnels, where it serves as the test medium, and in shock tubes, in which it is used as the 
driver gas. The perfect gas relations for helium at Mach numbers from 1 to 30 were 
presented in reference 1. Arave (ref. 2) presented the equilibrium flow properties for 
helium including the effect of first ionization for Mach numbers from 14 to 30 and quies- 
cent gas pressures of 10, 1, and 0.1 millimeters of mercury at a quiescent gas tempera- 
ture of 300° K. The quiescent gas is that gas, upstream of the incident shock, which is at 
rest with respect to the laboratory. With the advent of the expansion tube (ref. 3), it 
became necessary to know the equilibrium hypersonic flow properties of helium corre- 
sponding to quiescent gas pressures as low as 10“7 atmosphere. The present work was 
undertaken to satisfy this need and to include the effects of second ionization and particle 
interactions. 

Charts are presented of the thermodynamic properties of helium for temperatures 
to 100 000° K and pressures from 10“7 to 10^ atmospheres. The properties behind the 
incident, standing, and reflected shocks and at the in-flight and shock-tube stagnation 
conditions are also presented for incident-shock Mach numbers from 1 to 30, quiescent 
gas pressures from 10“7 to 1 atmosphere, and a quiescent gas temperature of 300° K. 
In-flight stagnation conditions result from an isentropic compression which brings the 
flow to rest, relative to the laboratory (or an object flying in an atmosphere), downstream 



of the incident shock. Shock-tube stagnation conditions result from an isentropic com- 
pression which brings the flow to rest, relative to the laboratory, downstream of the 
standing shock. 


SYMBOLS 

a speed of sound, centimeters/second 

a 0 radius of first Bohr orbit, 5.292 x 10 - 9 centimeter; speed of sound at p Q 

and T q , 97 242 centimeters/second 

B(i) electronic partition function of ith component 

Cp,Cy specific heats, at constant pressure and constant volume, respectively, 

ergs/mole -°K 

E energy per unit mass or per original mole of nonionized gas, ergs/gram 

E(i) energy per unit mass or per original mole of ith component, ergs/gram 

E n (i) energy of nth electronic level of ith component, °K 

e charge of electron, 4.80286 X 10”^ electrostatic unit 

g n (i) degeneracy of nth electronic level of ith component 

H enthalpy per unit mass or per original mole of nonionized gas, ergs/gram 

h Planck constant, 6.62517 x 10“^ erg-seconds 

I + ,I ++ ionization potential for neutral and singly ionized atoms, °K 

k Boltzmann constant, 1.38044 X 10“16 erg/°K 

M Mach number 

mjj mass of particle of kth component, m e = 9.1083 x 10“28 g ram anc j 

mjjg = m jjg+ = m He ++ = 6478 x 10 ^ gram 
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K p,i 

N(i) 


P 

p o 

^coul 


Q(i) 

R 

r 

A o,z 

s 

T 

To 

U r 

U s 

u 


u 2 


w. 


i» n 


equilibrium constant of jth reaction, dynes/centimeter^ 
number density of ith component, centimeter"^ 
integer indicating electronic level 

quantum number for electronic energy levels of ith component 
limiting quantum number for electronic energy levels of ith component 
pressure, dynes /centimeter^ 

reference pressure, 1 atmosphere (l. 01325 x 10® dynes/ centimeter^) 

pressure-correction term accounting for electrostatic forces, 
dynes j centimeter ^ 

partition function of ith component 

gas constant, 3.32824 x 10** ergs/mole-°K 

average distance from nucleus to nearest ion with charge ze, centimeters 

entropy per unit mass or per original mole of nonionized gas, ergs/mole-°K 

absolute temperature, °K 

reference temperature, 273.16 °K 

reflected-shock speed, centimeters/second 

incident -shock speed, centimeters/second 

gas velocity, centimeters/second 

gas velocity relative to incident shock, U g - U2, centimeters/second 
weighting factor 
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X(i) mole fraction of ith component 

Z compressibility factor 

z integer indicating charge in units of e 

a moles of He + per original mole of nonionized gas 

/3 moles of He ++ per original mole of nonionized gas 

y ratio of specific heats, Cp^Cy 

Ae + ,Ae ++ ionization-potential correction terms accounting for electrostatic 
interactions, °K 

p density, grams / centimeter^ 

p Q density at p Q and T Q , 1.7859 x 10"^ gram/centimeter ^ 

Subscripts: 

1 conditions ahead of incident shock (Tj = 300° K; pj 

2 conditions behind incident shock 

3 conditions behind standing shock in shock tube 

R conditions behind reflected shock 

sf in-flight stagnation conditions 

st shock-tube stagnation conditions 

ANALYSIS 

The thermodynamic properties of helium were determined from statistical 
mechanics in the same manner (with one relatively minor exception which is discussed 
subsequently) as reference 4. The reader is referred to this reference for all details. 
A compendium of the pertinent equations and spectroscopic constants (obtained from 


as indicated) 
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ref. 5) is presented in appendix A. These equations were programmed in FORTRAN IV 
language for solution on the IBM 7094 electronic data processing system. The computer 
program called, for temperature and pressure as input parameters and provided as output 
the following quantities: density, compressibility, mole fractions of individual species, 
energy, enthalpy, entropy, specific heats at constant volume and pressure, speed of sound, 
and the Debye-Hiickel pressure correction. 

The only difference between the properties presented in reference 4 and those pre- 
sented herein is a result of a slightly different treatment of the effect of particle inter- 
actions on the ionization potentials of the atom and singly charged ion. In the present 
work the lowering of the ionization potentials as predicted by Unsbld’s theory (refs. 6 
and 7) is given by 


and 


Ae + 


157 800 

M 2 ’ 


°K 


Ae' 


631 400 Orr 

K ) 2 ’ 


(1) 

( 2 ) 


for first and second ionization, respectively. The principal quantum numbers at which 
the combined potential of the nucleus and the nearest ion reaches a maximum were 
obtained from the equations 


1 3.33a 0 4.3a 0 

j 

and 

1 _ 2 - 4 4a 0 3a 0 

(n 2 *) 2 = r o,l +r o,2 


(3) 

(4) 


where a 0 is the radius of the first Bohr orbit (5.292 x 10"® cm) and r 0 z is the aver- 
age distance from the nucleus to the nearest ion with charge ze. The average distances 
were given approximately by Unsold 


and 



(5) 


cm 


r o 2 = 0.62 N(He ++ ) " 1//3 


( 6 ) 
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Here N(He + ) and N(He ++ ) are the average number of particles of the species He + and He ++ , 
respectively, per unit volume. Lick and Emmons in reference 4 approximated the fore- 
going set of formulas by the expressions 


Ae + = 5.65 x 10- 3 jN(e")] 7 , cm' 1 

and 

Ae ++ = 8.96 X 10" 3 [N(e')] 1//3 , cm -1 


( 7 ) 

(8) 


where N(e~) is the electron density. Because of the approximate formulation of the effect 
of particle interactions on the ionization potentials, there is little to recommend the 
choice of one calculation over the other. 

The program for the equilibrium thermodynamic properties of helium was com- 
bined with the program for calculation of normal -shock and stagnation-point conditions 
described in reference 8. With the two programs coupled, it was possible to compute the 
conditions behind an incident shock (a shock moving into a quiescent gas), a standing 
shock (a bow shock in front of a stationary body after passage of an incident shock), and 
a reflected shock (the reflection of an incident shock from a stationary, perpendicular 
surface) as well as the isentropic stagnation conditions upstream and downstream of a 
standing shock. 

A brief summary of the normal-shock equations is presented in appendix B. 


PRESENTATION OF RESULTS 


The thermodynamic properties of helium are presented as functions of temperature 
(to 100 000° K) and pressure (from 10 -r7 to 10^ atm) in figures 1 to 11. Values of the 
Debye-Hiickel pressure correction P cou i/P 0 for helium are given in the following table: 


T, 

°K 


Value of P coul /p 0 for - 


p/p 0 = 10- 2 

P/P 0 = 10° 

p/p 0 = i° 2 

p/p 0 = 10 4 

10 x 10 3 

5.08 X 10" 8 

5.11 X lO- 6 

5.16 x 10- 4 

5.27 X 10~ 2 

20 

3.12 X 10“ 6 

2.33 X 10-3 

3.95 x 10" 1 

4.77 X 10' 1 

30 

1.37 

1.35 

1.33 x 10° 

3.83 X 102 

40 

8.65 X 10- 7 

7.83 X 10" 4 

7.83 x 10" 1 

6.39 

50 

5.33 

5.60 

5.00 

5.14 

60 

3.70 

3.74 

3.77 

3.57 

70 

2.72 

2.72 

2.85 

2.69 

80 

2.08 

2.08 

2.13 

2.16 

90 

1.64 

1.64 

1.66 

1.74 

100 

1,33 

1.33 

1.34 

1.40 
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The values of the pressure-correction term are always small in the range of tempera- 
tures and pressures covered by the shock-tube calculations and, consequently, have been 
neglected in these calculations. 

The flow properties behind the incident, standing, and reflected shocks and at the 
in-flight and shock-tube stagnation points are presented in figures 12 to 28 for incident- 
shock Mach numbers from 1 to 30, quiescent gas pressures from lO"" 7 to 1 atmosphere, 
and a quiescent gas temperature of 300° K. 

All curves presented in the figures were plotted by a plotter which draws a straight 
line between succeeding data points.. 

SUMMARY OF RESULTS 

The thermodynamic properties presented herein are in good agreement with the 
less extensive results of reference 4. The largest differences, which are caused by the 
slightly different treatment of particle interactions, occur at the largest values of tem- 
perature and pressure considered in reference 4. For example, at a temperature of 
50 000° K and a pressure of 10^ atmospheres, the difference in the mole fraction of free 
electrons between the present report and reference 4 is about 0.003, or less than 
1 percent. 

The flow properties behind the incident, standing, and reflected shocks as well as at 
the in-flight and shock-tube stagnation points agree well with the more limited results 
presented in reference 2 except at the upper limit of the incident-shock Mach number 
range, where second ionization (not accounted for in ref. 2) becomes important. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 25, 1968, 

129-01-03-07-23. 
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APPENDIX A 


THERMODYNAMIC PROPERTIES OF HELIUM 

The equations used in the calculations for the thermodynamic properties are pre- 
sented in this appendix. 


Partition Functions 

The partition functions per particle for the individual species are 



Q(e-) = 2fe£j V2 l«^ 


where 


i 

B < He > ' I «'l,nSn< He > e 


-E n (He)/T 


n=l 

n 2 


B(He + ) - j w 2 g n (He+)e' En(He ^ 

n=l 


(A2) 


with 

nj = 8, if nj* £ 7 

nj equal to the largest integer less than (n-j* + lj, if n-|* < 7 
D-2 = 7, if n2* = 6 
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ri 2 equal to the largest integer less than ^12* + l), if n 2 * < 6 
Wj n = 1.0, if n S n*, for i equal to 1 or 2 


and 


w. 


i n = ni - nj*, if n = nj > n^, for i equal to 1 or 2 


Also, 


I + = 285 287 - Ae + 


I ++ = 631400 - Ae ++l 


(A3) 


The reductions in the ionization potentials are obtained from the equations 


where 


in which 


with 


A „+ 157 800 

Ae = -7 s- 

(” 0 2 

++ 631 400 




Ae 


K ) 2 


J 


1 3.33a Q 4.3a 0 




( n l) r o,l r o,2 


J 2 44a o , 3a n 


(n/) 2 r o,l r o,2 


J 


r Q) i = 0.62|N(He + )] 
r o,2 = 0.62[N(He ++ )] 


-1/3 


-1/3 


J 


N(He + ) = 2.687 x 10 19 Z^x(He + ) 
N(He ++ ) = 2.687 X 10 19 Z^x(He ++ ) 


(A4) 


(A5) 


(A6) 


(A7) 
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Equilibrium Composition 
The equilibrium constants for the two reactions 

He He + + e" 


and 


are, respectively, 


He + = He ++ + e 


K p,l _ Q(He + )Q(e~) a(a +2/3) 

P ” Q(He) “ (1 - a - / 3)(1 + a + 2/3) 

and 

K p,2 _ Q(He ++ )Q(e~) _ fifat + 20) 

P Q(He + ) Q!(l + a + 2/3) 

where 

Z = 1 + a + 2/3 
X(He) = 1 ~ ^ 

X(He+) = | l 

X(He ++ ) = | 

X(e-) = SL+2g 


(A8) 


(A9) 


(A10) 


With values of Ae + and Ae ++ given, equations (A8) and (A9) can be solved simultane- 
ously for a + 2(3 and the species mole fractions. Equations (A4) to (A7) can then be 
used to compute new values of Ae + and Ae ++ . The procedure is repeated with the new 
values until the conditions 


Ae k + 1 - Ae l 




10” 5 


Ktl - HI < 10 -5 

T 




(All) 


J 


are satisfied. Here the subscripts on Ae + and Ae ++ indicate the trial number. 
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Thermodynamic Properties of the Mixture 
The thermodynamic properties of the mixture are given by the equations 



and 

-S-=z 

RT 0 

where 


^ w 1 , n E n( He ) S n ( He ) e n 

3 t n=l 

2 TB(He) 


E(He) _ t 
RT 0 " T 0 


X(He)S|!l + x(He*)^} ♦ x(He~)=|^ + x(e-)fM 


E(He + ) _ T 3 
RT 0 T 0 2 


£ w 2 ,n E n(He + )g n (He + )e ^ 

n=l 


TB(He + ) 


D 

T 


E(He ++ ) = Tl|3 + 1+ + I ++ 
RT 0 T 0 [2 T 


E(e~) _ 3 T 
RT 0 2 T 0 



— = Tr> 

R ° 


)(h/rt q ) 


8T 


(A12) 

(A13) 

(A14) 

(A15) 

(A16) 

(A17) 

(A18) 

(A19) 

(A20) 
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and 


c v _c p t7 l + (T/z)(az/8 T ) p 
R R " Z 1 + (p/Z)(8Z/9p) T 


a r 


3yT\V 2 

5T 0 y 


1 



(A21) 

(A22) 

(A23) 


The derivatives (8/ 9T)p and (9/ 9p),p were evaluated numerically. The evalua- 
tion of these derivatives was necessary because of the unknown dependence of Ae + and 
Ae ++ on temperature and pressure. The Debye-Htickel pressure correction was obtained 
from the equation 


7.6591 X 
Po 


3/2/i \ 1/2 


"t) & 


(t) x ( e ") x ( He+ ) + 4X ( He++ ) 


(A24) 


The actual pressure is then given by 


p_ _ z (P\f T\ + p coul 
p o \ p oJVo) p o 


(A25) 


Because the correction was always found to be small, it was not included in the normal- 
shock calculations of this report. 


Spectroscopic Constants 

The spectroscopic constants for neutral and singly ionized helium are presented in 
the following tables. 


Neutral helium: 


Electronic level n 

g n ( He ) 

E n (He), 

°K 

1 

1 

0 

2 

3 

230025 

3 

13 

243 308 

4 

36 

267219 

5 

64 

275 354 

6 

64 

278975 

7 . . , 

64 

280 934 

8 

64 

! 282113 


Singly ionized helium: 


Electronic level n 

SnM 

E nM> 

°K 

1 

2 

0 

2 

8 

473 929 

3 

14 

561142 

4 

22 

592117 

5 

32 

606 633 

6 

44 

614 046 

7 

58 j 

618 700 
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NORMAL- SHOCK EQUATIONS 
Incident Shock 

A plane incident shock moving into a quiescent gas is shown in the following sketch: 



The conservation equations for the incident shock may be written as 

^2 = p l U s 

P 2 + p 2^2 2 = Pi + p l U s 2 

and 

H 2 + |3 2 2 . Hi + |u s 2 
Equations (Bl) to (B3) together with the equation of state 


(Bl) 

(B2) 

(B3) 


P 2 = p(p 2 ,H 2 ) 

form the system of equations for the unknowns P 2 , P 2 > and H 2 . 


(B4) 


Standing Shock 

A standing shock, such as a bow shock in front of a body, is shown in the following 
sketch: 
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The conservation equations across the normal portion of the standing shock are 


p 2 u 2 - P3U3 

p 2 + P 2 u 2 2 = p 3 + p 3 u 3 2 


and 


(B5) 

(B6) 


H 2 + |u 2 2 = H 3 + |u 3 2 (B7) 

The subscript 3 denotes quantities immediately downstream of the standing shock. The 
equation of state 

P3 = ^(p3’^3) (B8) 

completes the system. 


Reflected Shock 

A normal shock reflected from the closed end of a (one -dimensional) tube is shown 
in the following sketch: 


/ 
/ 
/ 
/ 
/ 


u r< 




u 2 




u R = ° 


/"Z7V ////// 7/7 


7/77-7-7 77 


The conservation equations when written in a shock-fixed coordinate system are 


^2 ( u 2 + U r) = PR V r 

(B9) 

P2 + P2^2 + U r) 2 = PR + PR U r 2 

(BIO) 

H 2 + |(u 2 + U r ) 2 = H R + |u r 2 

(Bll) 
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The equation of state is used in the form 

P R = P ( p R’ H r) 


(B12) 


Stagnation Conditions 

The stagnation conditions upstream of a standing shock (also called the "in-flight 
stagnation conditions") can be determined from the equation of state 


p sf = P( s sf> H sf) 


(B13) 


and the conditions 

Hsf = H 1 + | U s 2 " 

s sf = s 2 [ 


(B14) 


Similarly, the stagnation conditions downstream of a standing shock (also called 
the "shock-tube stagnation conditions") can be determined from the equation of state 


p st = p ( s st> H st) 

and the conditions 

H st = h 2 + | u 2 2 

> 

s st = s 3 


(B15) 


(B16) 


The numerical procedures required to solve the equations of this and the preceding 
sections of this appendix are fully detailed in reference 8. 
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Temperature, °K 


Figure 3.- Mole fraction of singly ionized helium. 
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Temperature, °K 


Figure 4.- Mole fraction of doubly ionized helium. 
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Temperature, °K 

Figure 10.- Specific heat at constant volume of helium. 
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